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Abstract. We report the discovery of a Lyman a emitting galaxy at z = 6.518. The single line was found in the 43 arcmin 2 
VLT/FORS field by slitless spectroscopy limited to the atmospheric window at A ~ 9100 A (sampling Lya- at 6.4 < z < 6.6). Its 
counterpart is undetected in a deep / band image and the line has an asymmetric appearance in a deeper follow-up spectrum. 
There are no plausible line identifications except for Lya with a flux of 1.9xl0~ 17 erg cnT 2 s and rest frame equivalent width 
> 80 A. The lower limit to the star formation rate density at z = 6.5 derived from our complete sample is 5x 10~ 4 M Q yr Mpc~ 3 , 
consistent with measurements in the Subaru Deep Field and Hubble Ultra Deep Field but approximately ten times higher than 
in the Large Area Lyman Alpha survey. This Lya emitter is among the very small sample of highest redshift galaxies known. 
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1. Introduction 

While reaching for the stars, our arms are still growing longer. 
Until a few yea rs ago, the mos t distant objects known were 
QSOs at z ~ 5 .8 ( Fan et al. 2000) and a handful of Lya emitters 
(LAEs) at 5 < z < 6 (Stern & Spinrad 1999, for a review). 
However, advances in the past years with respect to the red 
sensitivity of CCDs have enabled the detection of galaxies at 
even higher redshifts l lTaniguchi et alJ2 003. for a review). 

The detection of z > 6 galaxies allows to study the star for- 
mation rate in the early Universe, the modes of early galaxy 
formation and the interplay between the first galaxies and the 
intergalactic medium (IGM). With the large sample of Lyman 
break galaxies (LBGs) now available at z ~ 3 and z ~ 4, the 
star formation rate density (SFRD) at these redshifts is quite 
well determined (e.g. Steidel et al. 1999). At higher redshifts, 
there have been tentative determinations of the SFRD recently, 
but these are still based on poor statistics and some uncertain- 
ties in the m odeling ( H eavens et al. 2004; Bun ker et alJ l2004: 
iRicotti etaDl2004l) . At high redshift, we may encounter stars 
formed out of primordial elements, which differ strongly from 
those formed from enriched material, due to differences in ini- 
tial mass function (IMF) and stellar temperatures. These so- 
called Population III (PopIII) stars are responsible f or the pro- 
duction of the first m etals in the Universe llScannapieco et al.l 
l2003USchaere3l2003l) . The redshift range 6 < z < 7 is a very 
intriguing time during cosmic evolution, when hydrogen reion- 



* Based on observations obtained at the European Southern 
Observatory using the ESO Very Large Telescope on Cerro Paranal 
(ESO programs 066.A-0122, 068.A-0182, 272.A-5007). 



ization is believed to be basically complete and the IGM starts 
to be polluted with metals. Even though it is now believed 
that reionization o f most of the hydrogen takes place earlier 
(Kog ut et al.l l2003). large scale structure may inhibit a homo- 
geneous distribution of ionizing ra diation, introducing a large 
cosmic variance in the reionization (Gnedin 2000t lCiardi et alJ 
120031) . 

As the Lyo- line can be a prominent feature of young star 
forming galaxies, this line has been used extensively to search 
for distant galaxies. Blind spectroscopic searches are carried 
out from Earth ( e.g. lEllis et alfcOOll Martin & Sawickil2 004) 
and from space (Rhoads et al. 2004a; Pirz kal et al.l l2004). but 
most surveys are performed by narrow band imaging with 
ground based telescopes. Using narrow band filters sensitive 
to wavelengths in the atmospheric window at A ~ 0.9 pm, cor- 
responding to Lya emission at z ~ 6.5, and subsequent spec- 
troscopy, three groups have detected galaxies at z = 6.5. The 
first of these was found d ue to the lensing amplification (4.5x) 
by the cluster Abell 370 JHu et al.1120021) . two were found by 
surveying the huge field of view (814 arcmin 2 ) of the Subaru 
Suprime-Cam instrument (IKodaira et alJl2003l) and one during 
the course of the L arge Area Lyman Alpha survey (LALA, 
Rhoads et al. 2004b). Detections of even high er redshi f t, very 
strongly l ensed, galaxies ar e reported by Kne ib et alJ (|2004, 
Z = 7) and lPell6etal1d2004 z = 10). 

We have followed a different method to find Lya emitting 
galaxies at z ~ 6.5 and positively identify one galaxy at z — 
6.518. Throughout this Letter, we adopt a flat Universe with 
Q. M = 0.3, Q A = 0.7 and#o = 65 kms 1 Mpc -1 , 
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2. Observations 

Our approach to find LAEs at z ~ 6.5 is to use ground based 
slitless spectroscopy in combination with a medium band filter 
and imaging carried out with FORS2 at the VLT of a 7'x7' 
field with very low Galactic extinction (£Vi/=0.002) centered 
atRA = 04:01:55.9, Dec = -37:43:34 (J2000). The filter (ESO's 
z_SPECIAL+43, z sp from now on) has a central wavelength of 
9135 A and FWHM of 193 A and is therefore sensitive to Lya 
originating from z = 6.433 to z - 6.592. The filter reduces 
the sky background and overlap problems by providing short 
spectra. This technique provides a completeness limit on line 
flux and the redshift of lines with enough resolving power to 
separate possible contaminants from real LAEs. 

We used the holographic 600z grism which delivers a res- 
olution of about 1400 for 1" objects. Over a period of four 
months in the winter of 2002/03, an exposure of 7.6 hours was 
collected. Weather conditions were excellent resulting in 077 
seeing. Direct imaging was carried out through the z sp filter, to 
determine the sources of emission lines detected in the spec- 
troscopy, and Bessel-7 filter, to determine the / - z sp colour. A 
Gunn-v image was available from earlier observations of this 
field. Total exposure times were 6.33, 4.0 and 1.0 hours with 
seeing discs of 077, 077 and Of.' 6 for the z sp , I and v bands, 
respectively. 

Reduction was carried out in a standard way. Flux calibra- 
tion of the spectroscopic image was provided by observations 
of two spectrophotometric standards. For the z sp filter, the lit- 
erature values were convolved with the sensitivity of the filter 
to obtain the magnitude. Three sigma limiting AB magnitudes 
for a 1" aperture are 27.0 and 27.5, respectively for the z sp 
and / band. Special calibration were carried out to determine 
the distortion pattern of the FORS2 instrument using the mov- 
able slitlets of FORS2 in direct and spectroscopic modes using 
seven slit configurations to cover the entire field. These cali- 
bration images allowed us to determine the region on the direct 
image where the source of a spectrum or line on the spectro- 
scopic image can be located to within ~1 pixel or 0725. 

In november 2003, we have also carried out 2.1 hours 
of multi object spectroscopy in Director's Discretionary Time 
(DDT) with FORS2 using the 1028z grism, which delivers a 
resolution of about 2600 for the 1 " slits used in the mask. The 
eight frames were offset by 3" in an ABBA sequence. The ob- 
servations were reduced in a standard way, by subtraction of 
the overscan region and flat fielding with screen flats. Cosmic 
rays were rejected by combining the two pairs of frames with 
the same offset. Wavelength calibration and rectification of the 
two-dimensional spectra were carried out by measuring the po- 
sition of about 10 sky lines. The error in the wavelength cali- 
bration is dominated by fringing which influences the determi- 
nation of the center of the sky lines and limits the wavelength 
accuracy to 0.2 A. 



3. Results 

In the slitless spectroscopic image, three authors (JK, AC and 
SSA) have independently searched for emission lines by vi- 
sual inspection. Automatic extraction has been attempted but 
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Fig. 1. Part of the slitless image, including the targeted Lyor 
candidate and an [O in] emitter at z — 0.8 at the top right. 
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Fig. 2. Images in the v, / and z sp band of the region where pos- 
sible counterparts of the emission line can be located (between 
the dotted lines). The catalog numbers of the three possible 
counterparts are indicated. Note that KCS 1 166 is only detected 
in Zsp ■ 





E o.o- 



9120 9140 9160 9180 

Wavelength (A) 

Fig. 3. One and two dimensional slitless spectrum of the LAE 
candidate after smoothing by a 3 pixel boxcar along the dis- 
persion direction to increase the s/n without reducing the res- 
olution significantly. The dotted lines indicate the extraction 
aperture (4 pixels or 170). The dashed line indicates the signal 
to noise per pixel (after smoothing). 

fails to detect emission lines superposed on continua and re- 
sults in many fake detections. The visual inspection has a 
completeness level which corresponds to line fluxes of 1.0 
to 2.0xl0~ 17 erg crrT 2 s depending on the (unknown) wave- 
length of the line (values stated are for a line at the maximum 
of the filter transmission and half of it). Fig. ^shows a small 
part of the slitless image, including a LAE candidate which has 
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a signal close to the completeness limit. In the resulting com- 
bined list, there are 161 lines, each of which can originate from 
a counterpart within a region of about 50"xl" on the direct 
image. Taking all objects detected by SExtractor in the z sp im- 
age into account, the emission lines are related to 321 possible 
counterparts in the z sp image. The average number of possible 
counterparts per emission line is therefore two, which is con- 
sistent with one of these being the real counterpart as randomly 
placed regions with the above specified size cover on average 
one object in our z sp image. For each possible counterpart we 
have determined the associated line flux from the spectroscopic 
image (which depends on the position of the counterpart). 

Among the detected emission lines, 37 lines were part of an 
identifiable line pair, either close (e.g. [S II]/i6716,6731) or far 
apart (e.g. Hfi at 4861 A and [Olll]A4959). The corresponding 
87 counterparts have been removed from the sample of possi- 
ble LAEs. As the [OII]i3726,3729 doublet is only marginally 
resolved in our slitless spectra, we have not eliminated candi- 
date [Oil] emitters. From the remaining 234 objects, we have 
selected candidate Lyor emitters based on the singular and iso- 
lated appearance of the line, and the absence of emission in the 
/ band or strong / - z sp break of the possible counterpart. There 
is only one counterpart not detected in the / band. This ob- 
ject, KCS 1 166, has a line flux very close to the completeness 
limit (of 2.0xl(T 17 erg crrT 2 s ) and an/- z sp colour > 3.2 or 
> 1 .0 after subtraction of its line flux. The corresponding emis- 
sion line has two other possible counterparts, which are almost 
equally probable to cause the line (as they are close in posi- 
tion and brightness) but do not show a spectral break between 
the z sp and / band (see Fig.[2]and Table[0. The images of the 
three possible counterparts are rather faint for a reliable deter- 
mination of their sizes, but they are all consistent with being 
unresolved. The appearance of the emission line of this LAE 
candidate is asymmetric but it is difficult to quantify this as the 
s/n is low (see Fig. [3}. 

To determine the actual counterpart of this emission line 
and to cover a larger wavelength range with higher resolution 
and s/n, we have carried out multi object slit spectroscopy. Fig. 
[4] shows the two dimensional spectrum obtained with the slit 
mask, while Figs. |5] and [6] show the one dimensional spectra 
of the three possible counterparts. These spectra show that the 
line detected in the slitless image is from KCS 1 166. The other 
two objects are also line emitters: KCS 1 173 is identified as an 
[Oil] emitter at z = 1.13 and KCS 1147 as an [OIII] emitter 
at z — 0.73. A Gaussian fit to the emission line of KCS 1166 
results in acentral wavelength of 9141. OAorz = 6. 518 if iden- 
tified with Lyo-. KCS 1166 was also fit with a Gaussian pro- 
file truncated completely on the blue side (to simulate absorp- 
tion) whereafter it was (as all fitted profiles) convolved with a 
Gaussian instrumental profile. This fit (also shown in Fig. 
results in a FWHM of 280 km s 1 and redshift of 6.516 (for the 
unabsorbed line). Parameters derived from the Gaussian profile 
fits are presented in Table |2] 

Subtracting the line flux of KCS 1166 from the flux 
measured in the z sp image, we obtain a line-free contin- 
uum of 1.6+1. Ixl0~ 20 erg crrr 2 s A -1 or AB magnitude 
of 27.3. In case the observed line can be identified with 
Lyo-, no continuum flux is expected blueward of the line 



Table 1. Possible counterparts of the emission line 



Source 



Coordinates (J2000) 



/ 



KCS 1166 04 0143.48 -37 45 09.1 >27.4 >28.3 >3.2 
KCS 1173 04 0143.84 -37 45 08.9 25.2 25.0 0.3 
KCS 1147 04 0143.98 -37 45 09.1 "27.4 25.5 0.0 

Magnitude limits are lcr, errors are 0.1, except for a where it is 1.0 
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Fig. 5. One and two dimensional slit spectrum of object KCS 
1166. The dotted lines indicate the borders of the extraction 
aperture (4 pixels or 1"). The dashed line indicates the signal 
to noise per pixel. Shown as insets are two overlayed fits: a 
Gaussian and blue-absorbed Gaussian, both convolved with the 
Gaussian instrumental profile. 
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Fig. 6. Slit spectrum of object KCS 1 173 (left) and KCS 1 147 
(right). The dashed line indicates the background divided by 
50. 



3.2±2.1xlO- 20 ergcirT 2 s 



and the derived continuum redward of the line would be 
1 A -1 . The latter continuum flux im- 
plies an observed EW of 600+400 A for the line. As the con- 
tinuum detection is only lcr and we have conservatively used 
the high continuum flux, this value can be considered a lower 
limit. We have also measured the continuum level in the spec- 
trum of KCS 1 1 66 below (9005 < A < 9030 A) and above (9 1 65 
< A < 9244 A) the Lye line, resulting in values of -1.7+1.9 and 
2.1+2.5xlO~ 20 erg crrr 2 s A , respectively. The difference 
of 3.8+3.1 xl0~ 20 erg ctrr 2 s A -1 is consistent with the con- 
tinuum level redward of the line determined from the imaging. 
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Fig. 4. Complete two dimensional background subtracted spectrum of the slit with the sensitivity curve of the z sp filter indicated. 
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Table 2. Properties of the emission lines observed in the slit 
spectra 



Line 



Source 



Flux a 



FWHM" EW, 



i) 



Lyor ,11216 


KCS 1166 


6.518 


1.9+0.1 


160 


>80 


Oil ,13726 


KCS 1173 


1.129 


0.5+0.1 


110 


10 


011,13729 


KCS 1173 


1.129 


0.6+0.1 


110 


15 


H£ ,14861 


KCS 1147 


0.727 


0.4+0.2 


80 


20 


0111,14959 


KCS 1147 


0.727 


0.5±0.1 


<100 


25 


0111,15007 


KCS 1147 


0.727 


1.4±0.1 


<100 


65 


"10~ 17 ergcnT 


2 s b kms~' 
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4. Line asymmetry 

Due to the absorption of Lya photons by intervening neutral 
hydrogen, Lya lines are expected to have an asymmetric ap- 
pearance in constrast with non-resonant lines. A reliable iden- 
tification of an emission line with Lya depends therefore criti- 
cally on the assessment of the line profile. We have used various 
methods to quantify the asymmetry of the line of KCS 1 166. 

First, we have applied the three tests presented in 
Rho ads et alJ ( 120031) . A model-free test is provided by sym- 
metrizing the spectrum about an assumed line center and mea- 
suring how well the symmetrized line fits to the observed line. 
The x 2 value obtained by this method depends on the amount 
of continuum included, which is relatively symmetric around 
the line and therefore reduces the value of x 2 - It also depends 
on the error estimation. We have used the rms estimated from 
a source-free region in the slitless image above and below the 
actual spectrum. We have measured the x 2 value for a range of 
assumed line centers in steps of 0. 1 A. The lowest reduced x 2 
value measured was 1.35 for a conservative 20 A wavelength 
range including some continuum emission and 1 .46 for a range 
of 10 A which includes the line only. These values imply likeli- 
hoods of 1 .4% and 2.4%, respectively, that the line is consistent 
with its symmetrized profile. Two further tests attempt to mea- 
sure the actual asymmetry. For both, the wavelengths where the 
flux peaks and where it has 10% of its peak value on the blue 
and red side are determined. The wavelengths are respectively: 
A p = 9140.8A, A h = 9136.0A and A r = 9147.2 A. The wave- 
length ratio is then defined as ax - (A r - A p )/(A p - Ab) and the 

flux ratio as a/ — \ ' fxdA/ f " fxdA. The resulting values are 
ax - 1.33 + 0.39 and a/ = 1.12 + 0.07, where we have assumed 
an error equal to the binsize (0.8 A) for the wavelength and the 
measured sigma for the flux. These values are consistent with 
those measur ed for other high redshift Lya lines, according to 
lRhoadsetai]j2003l) . 

In addition to these three tests, we have also measured the 
skewness of the line for a range of wavelength windows. The 
skewness or third moment represents the degree of asymmetry 
of a distribution and is positive for a line if the values above the 
central wavelength are more spread out than those below it. The 
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Fig. 7. The skewness of the line for a range of window widths 
is represented as a histogram with one sigma errorbars. The 
dashed line represents the skewness of a Gaussian line profile, 
completely absorbed on the blue side and convolved with the 
instrumental profile. Regions where the uncertainty is due to 
either sampling or noise are indicated. 



result is presented in Fig. where the uncertainties for each 
window are also shown. These uncertainties are a combination 
of noise and sampling. For small window widths (<7 pixels), 
the uncertainties in wavelength sampling dominate, while for 
large widths (>12 pixels), the uncertainties in flux dominate. 
In between (7-12 pixels), a positive skewness is measured, 
consistent with a Gaussian emission line completely absorbed 
on the blue side and convolved with the instrumental profile. 

The results of all tests are consistent with the line being 
asymmetric, strongly suggesting that the line has excess emis- 
sion on its red side as expected for a high redshift Lye line. 



5. Discussion 

The absence of flux below the emission line at 9141 A, espe- 
cially in the deep / band image (see Tabled and the observed 
asymmetry support the identification of this line with Lya- at 
z - 6.518. The flux ratio between the line-free z sp and / bands 
of ~2 (or ~4 for continuum redward of the line only) is con- 
sistent with strong attenuation by the Lyman forest at high red- 
shift, but difficult to reconcile with a 4000 A or Balmer break. 
The identification with other lines than Lya is very unlikely. 
As the slit spectrum covers wavelengths up to 9262 A, we can 
exclude the identification with [O III] based on the absence of 
other lines in the spectrum. For the same reason, the identifi- 
cation with Ha is improbable as the accompanying [N II] lines 
would have to be ten times fainter and the Ha EWo would be 
in excess of 400 A. We would have resolved line doublets as 
[O II], Mg 11 and C iv (as shown by the spectrum of KCS 1 173, 
Fig. |6j and, if identified with [Oil], the line would have an 
EWo > 240 A, which is a factor of two higher than the largest 
EWo observed in low redshift samples of emission line galax- 
ies (e.g. Gallego et al. 1996). In addition, the measured / - z sp 
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colour, after line subtraction, is more consistent with a Lyman 
break than a Balmer jump. 

We can derive the SFR for the Lya emitter ba sed on the 
conve rsions for Ha and UV continuum given in Kennicutt 
(1998) and assuming a Lya/Ha ratio of 8.7 ( Brocklehurst 
Il97ll) . The Lya luminosity of the line is l.lxlO 43 erg s _1 , 
which results in a SFRL yff of 10.0±0.5 M Q yr , while the con- 
tinuum luminosity density of 6.6xl0~ 28 erg s _1 Hz -1 results in 
a SFRtjv > 9+6 M yr~' . Note that the SFRuv conversion used 
is calibrated only for the range between 1500 and 2800 A, while 
our UV continuum is measured near 1200 A, which can add ad- 
ditional uncertainties. 

The useful surface area of the survey is 6'. 18x6'93 or 42.8 
square arcminutes. We do not take into account the surface area 
covered by spectra of bright objects which may inhibit detec- 
tion of underlying emission lines as this is only a very small 
fraction (about 1%). The surface density of the complete sam- 
ple of LAEs (with fluxes > 2 x 10~ 17 erg cirT 2 s~') in our field 
is 0.023^^ arcmin~ 2 . Given the comoving volume of 18460 
Mpc 3 for our survey, we derive a lower limit to the SFRD at 
z = 6.5 of 5xl0 -4 M yr~' Mpc~ 3 (derived from the Lya flux). 
Using a much larger field of view and including LAEs with 
fluxes > 0.90 x IP' 17 erg cm' 2 s~'. iKodaira et al. (2003} ar- 
rive at a comparable SFRD of 4.8xl0~ 4 /z .65 M yr~' Mpc 3 . 
This SFRD is about ten times higher than that estimated for 
the complete sample of LAEs at z - 6.5 with fluxes > 2 x 
10~ 17 erg cirT 2 s -1 in the LALA survey (Rhoads et al. 2004b). 
It is interesting to note that the SFRD derived from our obser- 
vations is comparable with that derived from the abundance of 
i'-ban d drop-outs in the Hubble Ultra Deep Field (Bunker et al. 
2004), which low value poses a challenge for models suggest- 
ing that the bulk of star forming galaxies that reionized the uni- 
verse lie at redshifts just beyond z — 6. Larger samples of LAEs 
at z — 6.5 are needed to measure the clustering of these galax- 
ies and the resulting cosmic variance to obtain a better view of 
the SFRD at z > 6. 

One of the major challenges in observational cosmology 

is the detection of the first stars. These enigmatic objects 

may be responsible for the ionization of neutral hydrogen in 

the universe and should be formed from primordial matter. 

Although recent WMAP results seem to place reionization at 

I " II 1 

z > 10 (Kogutetal. 2003), the existence of metal-free stars 

at lower redshift depends strongly on the distribution of su- 
pernova products. The detection of large optical depth neutral 
hydrogen in front of z ~ 6 QSOs (e.g. Pentericci et al. 2002) 
illustrates that the ionization process might have occured very 
inhomogenously throughout the universe. IScannapieco et al.l 
(2003) predict that metal-free stars can be formed down to 
z < 4 and might be found in presently known LAEs (see also 
iFosburv et al.ll2003l) . The Lya EW of PopIII obje cts is pre- 
dicted to be very high (400 - 800 A, ISchaererl2003l) . although 
the observed EW can be diminished by H I absorption and de- 
pends therefore on the geometry of the ionized and neutral gas. 
The lower limit to the EW of KCS 1 166 (80 A) is not particu- 
larly high but we stress that the unabsorbed EWo can be signif- 
icantly higher. Table [5] shows the EWos for all known z - 6.5 
LAEs, as computed assuming that the observed continuum flux 
(or its lcr upper limit) is only present redward of the emission 



Table 3. Properties of known Lya emitters at z — 6.5 



Name 



Flux a EW! ReF 



HCM 6A 6.556 

SDFJ132415.7+273058 6.541 

SDFJ132418.3+271455 6.578 

LALAJ142442.24+353400.2 6.535 

KCS 1166 6.518 



0.65 25 (1) 

2.06 21 (2) 

1.13 44 (2) 

2.26 >85 c (3) 

1.90 >80 c (4) 



'10 17 erg cm 



b A c lp- lower limit d (l ) lHuetalJ J2002h . (2) 
, m lRhoads et"ail <2004hh . (4) This work 



II i I 

Rhoads et al. (2004b) and the 

galaxy described in this work have the highest EWos and are 
therefore the most likely LAEs at z — 6.5 to contain metal-free 
stars. Given the present observing capabilities and the uncer- 
tain Lya/He n flux ratio for PopIII objects ( Schaerer 2003), it 
would be just feasible to detect a strong He n line, a unique fea- 
ture of metal-free stellar populations. Further investigation is 
needed to determine the nature of the stellar population of this 
and other z - 6.5 galaxies. 
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